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The microscopic transport model UrQMD and a micro+macro hybrid model are used to calcu- 
late direct photon spectra from U-fU-coUisions at iSiab ~ 35 AGeV as will be measured by the 
CBM-coUaboration at FAIR. In the hybrid model, the intermediate high-density part of the nuclear 
interaction is described with ideal 3-|-l-dimensional hydrodynamics. Different equations of state of 
the matter created in the heavy-ion collisions are investigated and the resulting spectra of direct 
photons are predicted. The emission patterns of direct photons in space and time are discussed. 



I. INTRODUCTION 

In order to understand the phases of compact stel- 
lar objects and phase transitions in the early universe, 
a deep understanding of the phase diagram of strongly 
interacting matter is needed. The experimental tool to 
gain these informations is heavy-ion physics, where large 
nuclei are collided. In those collisions, nuclear matter 
is heated to temperatures that are believed to be suffi- 
cient to cause a phase transition from a hadron gas to a 
state where partonic degrees of freedom are relevant (the 
Quark-Gluon-Plasma QGP). Also at lower energies, but 
very high baryon densities novel phenomena like color su- 
per conducting phases [l| or even Quarkyonic Matter 
are expected in contrast to the well known hadron gas. 

From lattice gauge theory calculations, it is deduced 
that the transition between the hadronic and QGP phase 
is a cross-over [1, 0] , if the baryochcmical potential /ib 
is sufficiently small. The transition temperature is ex- 
pected to be around Tc ~ 170 MeV [1]. At high fiQ, 
first principle calculations of the phase structure are no 
longer possible, but symmetry arguments suggest that 
the phase transition line ends at high and zero tem- 
perature r = as a first order phase transition. If so, 
then a critical end point must exist, at which the cross- 
over turns into a first order phase transition. 

The search for this critical end point is the main mo- 
tivation for current and planned experimental heavy- 
ion programs of the SHINE-experiment at the Super 
Proton Synchrotron (CERN-SPS) H] and the CBM- 
experiment at the Facility for Antiproton and Ion Reser- 
ach (FAIR) 0. The current low-energy run at the Rel- 
ativistic Heavy Ion Collider (BNL-RHIC) also addresses 
this issue 

The CBM-coUaboration will measure uranium- 
uranium collisions in a fixed-target experiment at 
incident energies of up to i?iab = 35 AGcV. The 
proposed features for this detector include, besides 
subsystems to measure hadronic observables. electro- 
magnetic calorimetry in order to measure photons. 

Electromagnetic probes such as Icptons and photons 



have the advantage that they do not suffer from final 
state interactions and hence carry all information from 
their production vertex into the detector • Therefore, 
they also carry information from all stages of the heavy- 
ion collision. Hadronic probes, on the other hand, are 
mostly produced very late in the collision or rescatter in 
the late stages, so that they carry only indirect infor- 
mation from the early stages, such as flow patterns and 
fluctuations. 

Previous calculations of direct photons from transport 
theory include work with UrQMD by Dumitru et al. [ll| 
and Bauchle et al. [l2| and with HSD by Bratkovskaya 
et al. [l^ . Hydrodynamics has been used in many direct 
photon calculations, see e.g. [l^ - [2l| . 

Most photons produced in heavy-ion collisions, how- 
ever, come from hadronic decays, predominantly the de- 
cay of the 7r° and jy. Due to the long lifetime of both 
particles, they decay far outside the collision zone. In 
order to extract information about the fireball from pho- 
tons, the decay photon contribution has to be subtracted. 
The remaining photons, those that do not come from 
decays, are called direct photons. Experimental tech- 
niques for the extraction of direct photon yields from the 
inclusive measurements include the direct estimation of 
the background via invariant mass-analysis of the pho- 
tons [13, [2^ , the analysis of interference patterns (using 
a Hanburry Brown- Twiss analysis) [ll] and the extrapo- 
lation of the spectra of massive virtual photons (low- mass 
dileptons) to masslcss photons (25| . 

In this work, we investigate the infiuence of different 
evolutions of the bulk medium on the direct photon pro- 
duction in U-l-U-collisions at i?Lab — 35 AGeV using the 
model established in [l^]. We compare calculations with 
hadronic degrees of freedom to calculations with a first 
order phase transition to a QGP and calculations with a 
cross-over to chirally restored and deconfincd matter. In 
Section ini we briefly present the model and parameters 
used for these calculations, and in Section IIIII we show 
and analyse our calculations for direct photon emission. 



II. THE MODEL 
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The microscopic transport model Ultrarelativistic 
Quantum Molecular Dynamics (UrQMD) I26l428l pro- 
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vides the setting in which the present work is cal- 
culated. Calculations with UrQMD using standard 
options provides a baseline calculation with hadronic 
and string degrees of freedom throughout the nucleus- 
nucleus-interaction and with vacuum properties of all 
particles. 

Uranium has an eccentricity of e = 0.27. However, 
for the present study, the uranium nuclei are initialized 
spherically, which provides an intrinsic averaging over dif- 
ferent possible alignments of the projectile and target 
nucleus. We have checked that for untriggercd collisions 
(i.e. without selection on the alignment of the nuclei) the 
results are not sensitive to this choice. 

The hybrid option which is new in version 3.3 (29l - 
[3ll . l36j allows to substitute the high-density part of a 
heavy- ion collision with an ideal 3-f 1-dimensional hy- 
drodynamic calculation, where different assumptions on 
the matter present in the collision can be tested by 
varying the Equation of State. Non-equilibrium initial 
state interactions and final state scatterings in the dilute 
hadronic medium, as well as decays are calculated within 
the microscopic transport model. 

The transition from the particle-based transport de- 
scription to the density-based hydrodynamic description 
happens when the initial baryon-currents have decou- 
pled. We estimate this to be when the initial nuclei 
have completely passed through each other, which in the 
case of the system investigated here (U+U-collisions at 
E\3b = 35 AGeV) is at ^start = 3.7 fm after the very 
first scatterings. At this point, all particles that have 
interacted or are newly produced arc used to calculate 
energy-density, baryon number-density and momentum- 
densities, which are then taken as input for the hydro- 
dynamic calculation. During this transition, the system 
is forced into local thermal equilibrium, regardless of the 
state of actual equilibration of the matter. 

After the hydrodynamic evolution has proceeded, the 
Cooper-Frye-formalism [3^ is used to couple the hydro- 
dynamic calculation to the late transport stage. This 
transition happens individually for each transverse slice, 
characterized by the same position along the beam di- 
rection, when all cells in that slice have diluted to be- 
low a critical energy density, whose exact value depends 
on the Equation of State employed (see Table The 
hypersurface on which the transition happens is piece- 
wise isocronous and non-continuous. However, since af- 
ter that transition scatterings and decays are calculated 
with UrQMD, the effects of that should be negligible. 
For more details on the mappings, the reader is referred 
to Petersen et al. [3l|, [s!] . 



A. Equations of State 

Four different calculations are compared in the cur- 
rent work. Firstly, we compare direct photon spectra 
from hadronic systems, calculated in pure cascade model 
and in the cascade-|-hydrodynamic hybrid model with 



EoS 


tcrit 


HG-EoS 


5eo 


X-EoS 


7eo 


BM-EoS 


5eo 



TABLE I: The critical energy densities for the mapping from 
hydrodynamics to transport theory for the various Equations 
of state, eo = 146 MeV/fm'^ is the nuclear ground state energy 
density. 

Hadron-Gas EoS (HG-EoS) [H. The HG-EoS has the 
same degrees of freedom as the transport phase and 
therefore provides excellent means to investigate the ef- 
fects of instant thermalization at the transition from 
transport to hydrodynamic description and to inves- 
tigate the effects of the different kinetic descriptions. 
The hybrid model is also used with a Chiral EoS {x- 
EoS) [3a] featuring a rapid cross-over to chirally restored 
and deconfined matter, which is compared to the afore- 
mentioned hadronic descriptions and calculations with 
a MIT-Bag Model EoS (BM-EoS) with a strong first- 
order phase transition. Both x-EoS and BM-EoS have 
their phase transition at vanishing chemical potential at 
around Tc w 170 MeV. 



B. Intermediate stage in cascade calculations 

From the considerations above, it is clear that the hy- 
brid model consists of three stages, the pre-equilibrium 
early stage, the hydrodynamic intermediate stage and 
the dilute late stage. When comparing the stages be- 
tween the different variations of the model and among 
each other, we assign the same division of stages to the 
pure cascade calculations. Along these lines we define 
the early stage as < i < 3.7 fm, the intermediate stage 
as 3.7 < < < 11 fm and the late stage after t = 11 fm. 



C. Photon emission sources 

In the present work, direct photon emission is treated 
as a perturbation on top of the evolution of the fireball, 
so the underlying evolution of the hadronic medium re- 
mains unaltered by the calculation of direct photon spec- 
tra. This is justified, because the ratio of electromagnetic 
cross-sections (producing the photons) and strong cross- 
sections (governing the evolution of the bulk) is very 
small. 

The most important hadronic channels for the pro- 
duction of direct photons are tttt — >■ 7/3 and irp — >■ 771, 
which both are implemented in the transport, as well 
as in the hydrodynamic phase. The cross-sections for 
cascade-calculations are taken from Kapusta et al. (l3 |. 
while the rates used for the hydrodynamic description 
have been parametrized by Turbide et al. jlGj . Since no 
thermal partonic interactions arc modelled in UrQMD, 



3 




Px [GeVj 



FIG. 1: (Color Online) Overall direct photon spectra from 
U+U-coUisions. Calculations in pure cascade mode are shown 
as red crosses, hybrid calculations with HG-EoS (see text) 
are shown as solid blue line, x-EoS-hybrid calculations as or- 
ange dashed line and BM-EoS-hybrid calculations are shown 
as purple dotted line. 

emission from a QGP-mcdium is only taken into account 
in the hydrodynamic part of the model. Several minor 
hadronic channels are only implemented in one of the two 
models, such as strange channels (e.g. Ktt — > 7/v*) which 
are only present in the hydrodynamic calculations, and 
77-channels (e.g. TTT] — > 77r) which are only present in the 
transport calculations. Earlier investigations with this 
model have shown those channels to contribute about 
equally, but not significantly, to the overall spectrum of 
direct photons. The complete list of channels and a de- 
tailed explanation of the calculation is provided in [l2| . 

III. RESULTS 

We start with a comparison of the overall direct photon 
spectra calculated with the various Equations of State 
and with the transport-only approach. From Figure [U 
one can clearly see that the transport-only (crosses) and 
HG-EoS (solid line) calculations give very similar results. 
The BM-EoS and x-EoS calculations, on the other hand, 
yield significantly higher spectra than both hadronic cal- 
culations. The reason for this difference will be discussed 
below. 

Next, we investigate the contributions of the early, in- 
termediate and late stages to the overall spectra for all 
variations of the model. In Figure [2l we show that the in- 
termediate stage, which is calculated with hydrodynam- 
ics in the hybrid calculations, dominates the emission in 
all cases, although the excess of this stage is less signif- 
icant in the hadronic calculations of pure transport and 
HG-EoS hybrid calculations. 

Let us now take a look at the intermediate stages and 
investigate the contribution of hadronic and partonic di- 
rect photon emissions to the spectra from x-EoS and BM- 
EoS hybrid calculations. In Figure [3l we show the total 
contributions of the intermediate stages and the part of 
it that comes from hadronic sources (in case of the x~ 
and BM-EoS-calculations). We note that the hadronic 



Calculation 


Tsiope [MeV] 


A [GeV-2] 


xVd.o.f. 


Transport 


198.0±6.6 


2.09±0.74 


0.559 


HG-EoS 


203.5±8.0 


2.98±1.16 


0.532 


X-EoS 


214.8±6.1 


7.57±2.02 


0.249 


BM-EoS 


200.8±5.8 


15.39±4.45 


0.291 



TABLE II: Results of exponential fits to the spectra in the 
range < px < 3.5 GeV. The fit function is f{p±) = 



contributions are very similar in all cases, and that the 
excess observed in calculations with partonic degrees of 
freedom comes from the partonic part of the system. Fig- 
ure m shows the relative contribution of the QGP-phase 
in hybrid calculations with BM- and x-EoS. The con- 
tribution remains fairly constant at transverse momenta 
p± > 0.5 GeV, and is at about 85 % for the x-EoS and 
at about 65 % for the BM-EoS-calculations. 

What is eye-catching with the spectra is that despite 
the obvious differences in magnitude, the slopes of the 
spectra are very similar. Indeed, a closer anaylsis re- 
veals all inverse slope parameters to be about Tsiope = 
200 MeV, with only the X^EoS calculations being slightly 
higher (see Table [HI. 

Since emission of direct photons happens throughout 
the heavy-ion collision, it is prudent to look at the emis- 
sion times of photons. We do so in Figure El showing the 
average emission time of direct photons plotted against 
their transverse momentum for the two most abundant 
channels in cascade-only calculations and the time evo- 
lution of direct photon emission for various transverse 
momentum bins. 

We can see that the average emission time is nearly 
constant at about (t) « 7.6 fm over the whole p^-range. 
It is interesting to note that at intermediate p± direct 
photons from TTTr-scattcrings have a significantly lower 
average emission time (tTnr) ~ 4 fm than those coming 
from TTyO-coUisions. 

The time evolution of direct photon emission in the 
lower part of Figure [S] shows that at intermediate trans- 
verse momentum 1.5 < p± < 2.5 GeV the emission has 
two peaks, one after roughly 3 fm and the other one at 
8 fm. The average emission time in this momentum re- 
gion is {t{p± sa 2 GeV)) = 7.6 fm. At high transverse 
momentum, 2.5 < p± < 3.5 GeV, no significant peak at 
t > 5 fm is present, but the tail of the distribution is still 
large enough to have the average emission time as high 
as {t{p± « 3 GeV)) = 6.4 fm. 

At very low p± < 0.5 GeV, the emission is dominated 
by TTTT-scatterings, and is symmetrically centered around 
its average value of {t{p± w GeV)) = 8.8 fm, as is 
the 7rp-dominated photon emission at low transverse mo- 
mentum 0.5 < p±_ < 1.5 GeV; whose average time is 
{t{pj_ « 1 GeV)) = 7.6 fm. 

When looking at the relative contributions of the dif- 
ferent emission times iem to the direct photon yield as a 
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FIG. 2: (Color Online) The contributions of the early (green dash-dotted lines), intermediate (blue dotted lines) and late 
(red dashed lines) stages to the overall direct photon spectra (solid black lines) separately for all four variations of the model. 
Calculations in transport mode are shown in the top left panel, hybrid calculations with HG-EoS on the top right, X"EoS hybrid 
calculations are in the bottom left, and BM-EoS hybrid calculations are in the bottom right panel. 
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FIG. 3: (Color Online) Comparison of the contributions of 
the intermediate stages for all variations of the model. In 
addition, the hadronic contributions to the intermediate-stage 
emissions in x-Eos and BM-EoS are shown in thin lines. 



FIG. 4: (Color Online) The fraction of direct photons emit- 
ted from the QGP-phase in hybrid calculations with x-EoS 
(dashed orange line) and BM-EoS (purple dotted line) as a 
function of transverse photon momentum. 



function of direct photon transverse momentum in Fig- 
ure m we can assert the picture presented above. At 
low transverse momenta p± < 2 GeV, emission from 
5 < tern < 10 fm dominates and constitutes about 50 % of 



the overall contribution. At higher p± > 2.5 GeV, how- 
ever, emission from the very early stages of the collision 
starts to dominate. The contribution of photons emitted 
earlier than t^rn = 5 fm contribute more than 60 % to the 
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FIG. 7: Distribution of direct photons against the baryon 
number density pB at the emission point for cascade calcu- 
lations for all photons (red solid line) and for photons with 
p± > 2.5 GeV (blue dotted line, scaled by 5 ■ lO"). 



FIG. 5: (Color Online) Average emission time of direct pho- 
tons in cascade calculations as function of p± (top) and time 
evolution of the photon emission for various transverse mo- 
menta (bottom). The time evolution patterns from high 
transverse momenta have been scaled for readability. 
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FIG. 6: (Color Online) Fraction of photons that are emitted 
from < t < 5 fm (red dashed line), 5 < t < 10 fm (blue 
dotted line), 10 < t < 15 fm (green dash-dotted line), 15 < 
t < 20 fm and t > 20 fm (violett dotted line) as a function of 
transverse momentum. 



direct photon spectrum at p± > 3 GeV. 

The baryon number density at the point of emission 
of direct photons is investigated in Figure [71 Most 
photons come from coUisions below pB < 0.5po, with 
Pq = 0.16 fm~'^ being the nuclear ground state baryon 
number density. At high transverse momenta, the pho- 
ton emission is shifted from the low-density region to- 
wards high densities. We therefore assume it to be safe to 
employ vacuum cross-sections and vacuum thermal rates 
which do not take into account high-density effects on 
the p-spectral function. 



IV. SUMMARY 

Direct photon emission from U+U-collisions at i?iab = 
35 AGeV with 6 < 5 fm has been calculated and anal- 
ysed within a microscopic transport model and in a 
micro+macro hybrid model, in which the high-density 
phase of the transport model has been replaced by ideal 
3-|-l-dimensional hydrodynamic calculations. In the hy- 
brid framework, three different Equations of State gov- 
erning the behaviour of the matter in the hydrodynamic 
model have been compared. 

We find that the presence of partonic degrees of free- 
dom significantly enhance the amount of direct photons 
emitted from the medium, but not so much the slope of 
the direct photon spectra. At the energy under investiga- 
tion, we find that this enhancement comes directly from 
the partonic phase and is not due to a prolonged hadronic 
lifetime or hadronic contributions from the mixed phase, 
which is present in the BM-EoS hybrid calculations. 

The average emission time of direct photons is nearly 
constant at {t) w 7.6 fm over a broad range of transverse 
momenta in the case of cascade-only calculations, and 
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also the spectra from hybrid model calculations are dom- 
inated by emission from the intermediate stages. The 
average emission time of photons originating from tttt- 
scatterings drop to lower times at transverse momenta 
larger than p± > 1.5 GeV, but since at this p_L-range 
this process is subleading to the dominant Trp-channel, 
this does not affect the overall average emission time. 
Emission is found to come dominantly from areas with 
low baryon number density. However, a trigger on high 
p± will allow to select photons from the early stage. 

V. OUTLOOK 

The predictions for FAIR-energies shown in this work 
lay the foundations for future work analysing the ratio 
between direct photon emission and decay photons. Fur- 
ther work with this model will include the calculation of 



direct photon spectra at RHIC-energy and an analysis of 
the effect of changing the interface-parameters between 
the transport and hydrodynamic models on the direct 
photon emission. 
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